Spin-phonon coupling, q-dependence of spin excitations and high-Tc 
superconductivity from band models. 
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An understanding of spin excitations in cuprates is essential since the mechanism of high-Tc 
superconductivity might be hnked to spin fluctuations. Band calculations for long " 1-dimensional" 
unit cells of La2Cu04 show that the coupling between antiferromagnetic spin waves and phonons is 
larger for distortions of oxygens than for Cu or La. When this result is applied to a 2-dimensional, 
free-electron like band, it leads to a q-dependent spin excitation spectrum in good agreement with 
recent experiments. It is argued that important parameters for spin-phonon coupling, which comes 
out from the comparison between experiment and theory, are relevant for the mechanism of super- 
conductivity, and are large enough to explain a high Tc- 
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Many measurements have revealed the existence of 
complex superstructures and pseudogaps in high-Tc cop- 
per oxides. A picture of spatial, stripe-like regions with 
excitations of charge and spin modulations emerged quite 
early from neutron scattering experiments 1]. Angular 
resolved photoemission confirms that the bandstructure 
agree well with calculations, where one dispersive band 
gives rise to an ubiquitous barrel-like Fermi-surface (FS) 
(or a FS-"arc" when displayed in l/4th of the Brilloin 
zone) 0, Q- However, this arc is truncated and survives 
only in the diagonal direction {kx ~ ky) at low temper- 
ature, T, [j|. Phonons enter also into this complexity, 
as shown by the softening of some phonon branches at 
certain dopings [5,, ^, J} , by isotope effects on the pseudo- 
gap near the temperature T* , on Tc fS] , and on the gap 
structure itself [9[. Non-commensurate, g-dependent spin 
excitations appear as side-spots in neutron scattering at 
(0.5-9,0.5) and (0.5,0.5-9), where q varies linearly as func- 
tion of doping, X, up to a saturation at x ~ 0.12 [i0|. Re- 
cent measurements have established that the underlying 
spin excitations have a characteristic " hour-glass" shaped 
(g, w)-dispersion, even at non-superconducting composi- 
tions [ill, [l3, 13 • Here, it is shown that the character- 
istic spin dispersion can be understood in terms of spin- 
phonon coupling (SPC), and secondly it is concluded that 
phonons promote equal spin pairing and a large As/ for 
spin fluctuations, and therefore SPC is important for the 
mechanism of superconductivity. 

Previous ab-initio band calculations for supercells con- 
taining 1-dimensional (ID) phonon and spin-wave mod- 
ulations in the CuO bond direction ([1,0,0]), show large 
SPC within the CuO plane of these systems |lj|- This 
means that an antiferromagnetic (AFM) wave of the cor- 
rect wave length and the proper phase is stronger when 
it co-exists with the phonon [l^]. These results, in com- 
bination with 2D free-electron like bands, have been used 
for modeling of many normal state properties of the 
high-Tc materials [1^ ^^. Phonon softening, dynam- 
ical stripes, correlation between q and x, smearing of the 
non-diagonal part of the FS, and abrupt disappearance of 
the spin fluctuations at a certain T*, are possible conse- 



qucncies of SPC within a rather conventional band. The 
present calculations consider SPC between spin waves 
and four different types of phonon distortions. 

First, we discuss the ab-initio band calculations made 
for La(2_2;)Baa;Cu04 (LBCO), where the virtual crys- 
tal approximation (VCA) is applied to La-sites to ac- 
count for the doping. The calculations are made for 
long supercells oriented along the CuO bond direction, 
by using the linear Muffin-Tin Orbital method (LMTO) 
in the local spin-density approximation (LSD A), as has 
been described previously [l5|. Phonon distortion am- 
plitudes (u) and the size of Cu moments (m) in spin 
waves are necessary input to these calculations. The T- 
dependences u^ « iksT/Ku and m^ « ksT/Km are 
vahd for not too low T. Here Kp — cPEldjp^, E is the 
total energy, and p — u or m, respectively [15|. The force 
contants K^ for the different atoms are taken from mea- 
surements on YBa2Cu307 [l3,[iM|- The resulting u/ao, 
where ag is the lattice constant, are shown in the Table 
for T « lOOi^ (the interesting temperature range for Tc 
and T*). An approximate calculation of Km for a short 
wave in HgBa2Cu04 corresponds to magnetic moments 
m on Cu of about r^ 0.09/iB at 100 K 19]. Such moments 
are typically obtained by application of a magnetic field 
of the order of ±5 mRy [155. 

Results of calculations are presented here showing a 
varying degree of SCP for different phonons. These cal- 
culations consider distortions in a cell with 16 formula 
units (8ao along x) with displacements of La and api- 
cal O along z, and Cu and planar-0 along S, and the 
phonons may co-exist with AFM spin-waves. For the 
latter phonon there is a positive SPC when the nodes of 
the AFM waves are located at the " compressed" Cu sites 
(when the 0-atoms are moved towards the Cu) . Optimal 
SPC for displacements of out-of-plane atoms (La and api- 
cal O) occurs when these atoms move towards the CuO 
plane near the region of the AFM node. The calculated 
results for the maximal potential shifts on Cu caused by 
phonons (V^) and by spin waves (VI" in the spinpolar- 
ized potential) for the 4 types of movements in the cell are 
shown in the Table, together with the result for V^ with- 



out phonon. That the phase between the phonon and 
spin wave is crucial for having a positive SPC is shown 
by the fact that V™ for an " apical-oxygen" phonon is re- 
duced from 10 to below 8 mRy when the phase between 
the phonon and the spin wave is shifted by tt. 

Measured [13, [ISl and calculated [2l|, [24I phonon fre- 
quencies for zone boundary phonons for YBa2Cu307 are 
indicative for typical hu) for each atomic character of the 
phonon DOS. Phonon energies and partial phonon DOS 
for each site, Ngite have been calculated by Chen and 
Callaway [2^ for a similar cuprate system, Nd2Cu04, 
and the same NgUe can be assumed to be representa- 
tive for LBCO. From these references it is seen that the 
main La modes are at 10-20 meV, Cu at 20-30 meV, pla- 
nar O near 50±2Q meV, and apical O at 60±15 meV 
[13, mi Hi. Our estimations of Nsite from ref. [H are 
shown in the last part of Table I. A rescaling is made 
for the wavelength. The LMTO results are calculated for 
short waves, 8ao, but the waves are longer at the doping 
X = 0.16 (corresponding to g w 0.08), when V^ and V™ 
are larger by factors of ~ 1.1 and ^1.5, respectively [Tg. 
These factors might be fine tuned later if the calculated 
wave lengths fall outside the expected range, but in the 

from. 



first calculations we calculate the total V* 



K' 



= ^(i.i*^,p, + i.5*y;,'j)*iv. 



(1) 



where i is the site index. The resulting V* are 17, 18, 23 
and 22 mRy at the energies centered around 15 (La), 25 
(Cu), 50 (plane-0) and 60 meV (apical-O), respectively. 
These values are used in the subsequent calculations of 
the q — Lu dispersion. As discussed later, spin waves at 
higher energy are independent of the phonons, and V^ = 
0. 

In a second step we use the parameters V* in a 2D 
nearly free-electron (NFE) model. The AFM spin ar- 
rangement on neighboring Cu along [1,0,0] in undoped 
LBCO corresponds to a potential perturbation, V{x) = 



TABLE I: Three first lines: Phonon distortion amplitudes, 
u/ao, induced potential shifts, Vq (mRy), and exchange split- 
ting, Vq" (mRy), caused by spinwaves {hbH — ±5 mRy). All 
shifts/splittings refer to the maximal value on a Cu-site, and 
they are determined in self-consistent LMTO calculations for 
supercells of length 8ao. Remaining lines: The partial charac- 
ters of the phonon DOS, N;, for different sites, i, as estimated 
from ref [2311. 



wave 


no-phon 


PLO:, 


La^ 


ap-O^ 


Cu, 


u/ao 


- 


0.014 


0.021 


0.017 


0.024 


Kf 


- 


15 


2 


5 


3.5 


KT 


8 


12.5 


12 


10 


8.5 


Npi-o 


- 


0.6 


0.0 


0.2 


0.2 


NLa 


- 


0.0 


0.65 


0.0 


0.35 


Nap_0 


- 


0.35 


0.0 


0.65 


0.0 


Ncn 


- 


0.2 


0.2 


0.0 


0.6 



V^exp{—iQ ■ x) (and equivalently with y along [0,1,0]). 
The periodicity in real space is defined through the Cu- 
Cu distance, and a gap of size 2V1* appears at the zone 
boundary, at Q/2. A further modulation (q) of this order 
into ID-stripes perpendicular to x (or "checkerboards" in 
2D along x and y) is achieved by a multiplication of the 
potential by exp{iq ■ x) , where q < Q. Totally, this makes 
V{x) = V^exp{—iQx ■ x), where Qx = Q — <!■ The peri- 
odicity of this potential is now larger, and the gap moves 
away from Q/2 to (Q — q)/2. Magnetic side spots appear 
at q/2 surrounding Q/2 in probes which can separate the 
two periodicities. A 3x3 eigenvalue problem with matrix 
elements Hii=E-kl- k^, H22 = E - [k^ - Qx? - fey, 
H33 - E~kl-{ky-QyY, H12 = i/i3 = V^ and ifsa = 0, 
is solved. The lowest FE band (when V_* = 0) contains 2 
electrons up to Ep ^ 0.15i?y, when the effective mass is 
one. The total band width is quite close to the real band. 
The model is entirely 2-D, (the DOS is constructed from 
a sum over all states in the plane), which is reasonable 
in view of the very small band dispersion along kz of real 
band structures for the cuprates fl4]. 

An important result of the 2D-NFE model is that it 
leads to a correlation between doping and the amplitude 
of V/i [la]. The reason is that the gap opens along (fej,, 0) 
and (0, fcy), but not in the diagonal direction. The com- 
bined effect is that the dip in the total DOS (at which 
Ep should fall for optimal doping) will not appear at the 
same band filling for a small and a wide gap, even if the 
q- vector is the same. Alternatively, since the gap should 
appear at the same energy (at Ep) for different VI*, one 
has to vary the q-vector until it fits. The present calcu- 
lations are made for 0.16 holes per Cu (which is close to 
x in ref. [l^) with V* from eq. 1. 

The results are shown in fig. 1, with positive and nega- 
tive q displayed symetrically, and the spectrum is shaped 
like an hour-glass with a "waist" at intermediate energy. 
The points below 70 meV are for the coupling to the 4 
types of phonons. The general shape of the q — to de- 
pendence and the amplitudes of q are similar as in ref. 
[la ], with the smallest q at about 0.05 for SPC domi- 
nated by 0-p phonons around 50 meV. The weaker SPC 
for La at low energy makes q larger. The calculated 
results depend on the different parameters in the NFE 
model. Larger band mass leads to smaller amplitudes of 
q. The variations of g is a result of the differences in 
V* at different energies. Consequently, if no mixing of 
the phonon modes (through the A^siie-coefficients) were 
made it would increase the variations of q, even though 
the general hour-glass shape would remain. 

The uppermost part of the spectrum needs some com- 
ments. Spin waves and phonons are tied together at 
the same frequency and g-vector for the optimal mecha- 
nism of SPC. Spin waves with higher frequency than the 
phonons cannot profit from this mechanism. Thus, it is 
assumed that spin waves with frequencies higher than the 
phonon frequencies are independent of lattice vibrations, 
whereby V^ = 0. This "phonon independent" result is 
put rather arbitrarily at huj ^^120 meV in fig. 1, which 



is about twice the highest phonon frequency. The cal- 
culated wavelengths for the SPC modes below '^70 meV 
are found to be longer than the length of the supercell in 
the LMTO calculation, which justifies the factors 1.1 and 
1.5 in eq. 1. However, the solution without SPC leads to 
a wave which is even shorter than the cell in the LMTO 
calculation. The the factor 1.5, in eq. 1, is not justified 
for this (phonon free) mode, and a corrected calculation 
with y* =8 mRy is more appropriate, with a solution for 
q = 0.24, see fig. 1. In contrast to the cases with SPC 
and larger VI*, this solution implies a very short wave, 
4-5 ao only, where V* might be reduced even more. This 
will, because of the self-consistent feed-back between spin 
density and potential, lead towards a vanishing spin wave 
|15j . Therefore, it is difficult to imagine larger q for spin 
excitations with small V* at this doping. 

Other high-energy solutions with small V* exist for 
non-equal q- vectors along x and y [la|, where the upper 
of two gaps corresponds to a doping of 0.16. These so- 
lutions are found within some range of q, but the dips 
in the DOS are relatively weak. A satisfactory solution 
is found for qx and qy near 0.14 and 0.11, respectively. 
The average of the two vectors, 0.125, is comparable to 
the wave length in the LMTO calculation (where V'"=8 
mRy). Two gaps remain in the DOS for larger separa- 
tion of the two g-vectors, but the dips are found at too 
high and too low doping. The existence of this multitude 
of solutions indicates that large broadening and damp- 
ing is expected at high energies. Such fluctuations are 
spread in energy and momentum. They decay rapidly, 
since they are not linked to the phonon spectrum. 

Less doping will make the waves longer and the waist 
in the q — huj diagram becomes narrower. The narrowing 
should be noticed at all energies even at the highest E 
where spinwaves are excited without the help of phonons. 
Smaller q-vectors have been observed recently in lightly 
doped Lai.96Sro.o4Cu04 [2J], in agreement with this pre- 
diction. However, the spin modulation has turned from 
parallel (the Cu-0 bond) to diagonal direction at such 
low doping. In addition, the waist in the observed spec- 
trum is found at lower energy than for parallel doping 
[24|, which in the SPC-model implies increased SPC for 
phonons at lower energy. Ab-initio calculations for mod- 
ulations along the diagonal direction are needed to verify 
if La- or Cu-modes grow in importance relative to the 
high energy 0-modes. 

Heavier 0-isotopes will decrease the frequencies for the 
phonons and the coupled spin waves, and move the waist 
to lower E. They may also decrease u, at least at low 
temperature, which in the model for SPC will make V* 
smaller. From the correlation between V?" and wave 
length it is expected that q becomes larger. In all, this 
leads to a wider waist at lower energy, while the up- 
per part of the spectrum should be insensitive to isotope 
shifts. 

Undoped cuprates are generally stable AFM insula- 
tors. The exchange enhancement for spin waves in doped 
cuprates is moderately large in the absense of phonons, 



Doping = 0.16 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

q (r.l.u.) 



FIG. 1: Calculated q — hto relation from the 2D-NFE model 
and the parameters Vq for doping x =0.16 as in the experi- 
ment of VignoUe et al [1^. The lines connecting the points 
indicate the expected continuation for intermediate energies. 
The solution without SPC, the uppermost point, is not very 
precise (see text). 



and it becomes peaked if a phonon can intervene. This 
can be translated into a large phonon deformation po- 
tential. The unusual part of it is that this deformation 
potential resides in the equal spin part, and hence it 
is more correct to associate this deformation potential 
with a large As/ for spin fluctuations. The calculated Vq- 
parameters determine the spin excitation spectrum for 
large SPC, but they are also the key parameters for es- 
timates of superconducting T^- For instance, V^ corre- 
sponds to the monopolar matrix element for electron- 
phonon coupling, Xph- The good agreement between the 
calculated spin excitations and experiment suggests that 
the calculated V^'s are of the correct order. It is not 
clear how the contribution to T^ from coexisting phonon 
and spin fluctuations will be shared, but a simple ac- 
count for separate Xph and As/ gives rather large val- 
ues for Tc- The estimates are based on A = NP /K^, 
where N is the DOS at Ep (« 0.5{eV -Cu- spiny^), and 
/ = (AVP/Au). From the Table and with oq = 3. 781, 
AV^/Au ~ 3.8eV/A for a pure (plane) 0-phonon, which 
makes Xph ~ 0.5, for K^ of the order 15 eV/A^. If 
phonons stir up a spin fluctuation (through SPC) one 
should only consider the equal spin part in AV"^, which 
is almost of the same order as AV^ for the 0-phonon (see 
Table). Moreover, one should add a Km to K^ because 
of the energy of the spin wave, but this contribution can 
be negative if the coexistence makes the phonon softer 
{Km without phonon may be large and positive). By as- 
suming Km ^ this gives As/ ~ 0.4, and it will increase 
faster than Xph when the doping becomes smaller. The 
BCS formula Tc = 1.13 • Huj ■ exp{—l/X) makes Tc of the 
order 80-40 K, for these two A's, when hui — 50'meV. 



This is very approximate, but it shows that the couphng 
strength can be sufficient for a large Tc- Note also that 
since the exchange enhancement in LSDA is too weak to 
stabilize AFM in undoped LSCO, it will also produce a 
too low value of As/ . 

The total SPC becomes too large at underdoping when 
the wave modulations are longest, so N (and A) are re- 
duced because of the pseudogap, and hence the pseudo- 
gap competes with superconductivity. On the overdoped 
side A is reduced because of a small AV, and Tc -^ 0. 
Optimal conditions for a high Tc are found at interme- 
diate doping, but the estimates are not precise enough 
for finding the exact value or the exact doping limits for 
vanishing Tc- The role of phonons decreases as a; ^ 0, 
since the ratio V^/V^ is growing at small doping. 

In conclusion, these calculations suggest that the hour- 
glass shaped dispersion of the spin wave spectrum is 
a consequence of different degrees of SPC for different 



phonon modes. The " half- breathing" 0-mode in the mid- 
part of the phonon spectrum is found to be the most effi- 
cient one, which explains the waist in the spin wave dis- 
persion. The SPC-NFE model simulates several normal 
state properties surprisingly well [16|. There is nothing 
special with the electronic band structure in this scenario. 
The LDA bands are essentially correct for doped systems, 
although the exchange enhancement is underestimated in 
normal LSDA calculations |25| . Finally, it is argued that 
the observed spin excitations provide important informa- 
tion about the mechanism of superconductivity, since the 
excitation spectrum can be related directly to the cou- 
pling parameters. In particular, there is an unusual con- 
tribution to the equal spin pairing parameter X^f coming 
from the coupling to phonons. 

I am grateful to B. Barbiellini and C. Berthod for var- 
ious discussions. 
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